ASYMPTOTIC BEHAVIOUR OF A RAPIDLY ROTATING FLUID WITH RANDOM 

STATIONARY SURFACE STRESS 



ANNE-LAURE DALIBARD* 

Abstract. The goal of this paper is to describe in mathematical terms the effect on the ocean circulation of a random 
stationary wind stress at the surface of the ocean. In order to avoid singular behaviour, non-resonance hypotheses are 
introduced, which ensure that the time frequencies of the wind-stress are different from that of the Earth rotation. We 
prove a convergence result for a three-dimensional Navier-Stokes-Coriolis system in a bounded domain, in the asymptotic 
of fast rotation and vanishing vertical viscosity, and we exhibit some random and stationary boundary layer profiles. At 
' last, an average equation is derived for the limit system in the case of the non-resonant torus. 
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1. Introduction. The goal of this paper is to study mathematically a problem arising in ocean 
dynamics, namely the behaviour of ocean currents under stimulation by the wind. Following the books 
CN , by Pedlosky [jQlEl] and Gill [9], the velocity of the fluid in the ocean, denoted by u, is described by the 

incompressible Navier-Stokes equations in three dimensions, in rotating coordinates, with Coriolis force 

CL 

■ p{dtu -f u • Vu + 2ne A u) - AhAhU - A^d^u = Vp, t> 0, (x, y, z) e U(t) C 

divu = 0. 

-4— > . 

, In the above equation. Ah and Az are respectively the horizontal and vertical turbulent viscosities, p is 

the pressure inside the fluid, p is the homogeneous and constant density, and fie is the rotation vector 
of the Earth > and e is a unitary vector, parallel to the pole axis, oriented from South to North). 
U{t) is an open set in R"^, depending on the time variable t: indeed, the interface between the ocean and 
the atmosphere may be moving, and is described in general by a free surface z = h{t). 
\^ ' In order to focus on the influence of the wind, let us now make a series of crude modehng hypotheses 

Q^ , on the boundary conditions: first, we assume that the lateral boundaries of the ocean are flat, and that 

■ the velocity u satisfies periodic boundary conditions in the horizontal variable. We also neglect the 

fluctuations of the free surface, namely, we assume that h{t) = aD, with a,D positive constants. This 
CO ' approximation, although highly unrealistic, is justifled by the fact that the behaviour of the fluid around 

, the surface is in general very turbulent. Hence, as emphasized in [6], only a modelization is tractable 

and meaningful. Let us also mention that the justification of this rigid lid approximation starting from 
a free surface is mainly open from a mathematical point of view: we refer to [T| for the derivation of 
Navier-type wall laws for the Laplace equation, under general assumptions on the interface, and to [14] 
, for some elements of justification in the case of the great lake equations. At last, we assume that the 

^ ' bottom of the ocean is fiat; the case of a nonfiat bottom has already been investigated by several authors, 

5^ , and we refer to [61, r8l,fT7j for more details regarding that point. 

As a consequence, we assume that U{t) = [0,aiH) x [0, a2-ff ) x [0, aD), where if > is the typical 
horizontal lengthscale, and u satisfies the following boundary conditions 

u is periodic in the horizontal variable with period [0,aiH) x [0, a2i?), 
W|z=o = (no slip condition at the bottom of the ocean), 
dzUh\z=aD = ^oo' (infiuence of the wind), 
U3,\z=aD = (no fiux condition at the surface). 

Let us now reduce the problem by scaling arguments. First, we neglect the effect of the horizontal 
component of the rotation vector e, which is classical in a geophysical framework (see [4]). Furthermore, 
we assume that the motion occurs at midlatitudes (far from the equator) , and on a "small" geographical 
zone, meaning H <^ Rq, where i?o is the earth radius. In this setting, it is legitimate to use the so-called 
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f -plane approximation (see [7]), and to neglect the fluctuations of the quantity 63 • e with respect to the 
latitude. In rescaled variables, the equation becomes 



dtu''-" + u'^'' ■ Vu''-' + ies A u^-"^ - J/^A^u^'"' - i^.d^X^'' + (^gf^) = 0, (1.1) 



£ 

where 



U ^ _Ah_ ^ LA, ^ D 



and U is the typical horizontal relative velocity of the fluid. We are interested in the limit 

j/^ < 1, e < 1, Vh ^ 1. 

Such a scaHng of parameters seems convenient for instance for the mesoscale eddies that have been 
observed in western Atlantic (see [20]). One has indeed 

C/-5cm-s"\ 100km, Z) - 4km and - 10" 

which leads to e ~ 5 x 10"^. Possible values for the turbulent viscosities given in [20j are 



Ah ^ 10^ kg - m ^ • s ^ and ^2 1 kg • ni 



s 



so that i^z ~ 10~^, i^/i ~ 1. In the rest of the article, we denote by f the small parameter Vz, and we 
assume that t'/i = 1. Additionally, we do not take into account the shallow water effect, and thus we take 
7 = 1, even though this is not consistent with the values of H and D given above. Indeed, the thin layer 
effect, which corresponds to 7 ^ 1, is expected to substantially complicate the analysis, but without 
modifying the definition of boundary layers. Thus, in order to focus on the infiuence of a random forcing, 
we study the classical rotating fiuids equation (see for instance [1]), that is 

dtu'^" + u'^" ■ Vu'-'' + iea A u"^" - /^hu'^" - vdlu'^" + Vp = 0. (1.2) 

e 

Moreover, the amplitude of the wind stress at the surface of the ocean may be very large; thus we set 

AoSoD 



U 



where Sq is the amplitude of the wind velocity, and we study the limit (3 ^ 00. Equation (|1.2p is thus 
supplemented with the boundary conditions 

9^<ta=/?^^ (1-3) 

tt„ = 0. 

A\z—a 

Additionally, u^''^ is assumed to be T^-periodic in the horizontal variable x/j, where M^/[0, ai) x 

[0, a2). In the rest of the paper, we set T := x (0, a). The assumptions on the wind-stress will be 
made clear later on. 

1.1. General results on rotating fluids. Let us now explain heuristically what is the expected 
form of u^'" at the limit. Assume for instance that v — e and that the family u^-" behaves in L?{[Q, T] x T) 
like some function u°{t,t/e,x), with u° € L°°{[0,T] x [0,oo) x T) sufficiently smooth. Then it is natural 
to expect that uP satisfies 

drU" + 63 A li° = 0, 

divu" = 0, (1.4) 



"3|z=0 ~ "3|z=a — 0- 
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In fact, the above equation can be derived rigorously from l|1.2p if the dependence of the function 
with respect to the fast time variable r is known a priori] the goal of the two-scale convergence theory, 
formaHzed by Gregoire Allaire in after an idea of Gabriel N'Guetseng (see [18]), is precisely to justify 
such derivations in the context of periodic functions. However, in this paper, we do not need to resort 
to such techniques; our aim is merely to build an approximate solution thanks to formal computations. 

In view of (|1.4p . we introduce the vector space 



H := {u e L^(T)^ divu = 0, 



-a = 0} 



We denote by P the orthogonal projection on TL in L^(T)'^, and we set L := P(e3 A-). Notice that P differs 
from the Leray projector in general, because of the no-flux conditions at the bottom and the surface of 
the fluid. It is known (see for instance jij) that there exists a hilbertian basis [Nk)ke-L^\{o} of H. such 
that for all k, 



(ea A Nk) = iXkNk with 



W\ 



where k' — (27rfci/ai, 271^2/02, tt/cs /a). The vector N)~ is given by 

Nk{xh,z) 



coslk'^z)ni(ky 
cos{k'^z)n2{k) 
sin(fc3z)n3(fc)^ 



where 



ni(fc) 
7^2 (fc) 
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y/aia2a\kf^\ 

—^—{-tk[ + k'^Xk) 
^/aTa^\kf,\ 



if kh + 0, 



^Ja\a2a\k'\ 



and 



n\{k) 
n2{k) = 



sgn(fc3) 



else. 



Let £(t) : H ^ H he the semi-group associated with equation (|1.4p . i.e. £(r) — exp(— tL) for r > 0. 
We infer from equation (|1.4p that u'^{t,T) G H almost everywhere, and that there exists a function 
such that 



Consequently the main effect of the Coriolis operator L is to create waves, propagating at frequencies 
of order e~^. The goal is now to identify the function u° j which in general depends on the slow time 
variable t. This is achieved thanks to filtering methods, which were introduced by S. Schochet in |23j, 
and further developed by E. Grenier in [10| in the context of rotating fluids. Precisely, setting 

= exp u^'" 



a that u^'"^ converges 



it is proved in [H [TT] in the case of Dirichlet boundary conditions at z = and z 

strongly in ^[^^^([0, 00) x T) towards a function u\. Moreover, the function u° satisfies a nonhnear 
equation of the type 



dtul + Q{utul) 



(1.5) 
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where the quadratic term (3(u°,u°) corresponds to the filtering of oscillations in the non-linear term 
u*^'"- Vm^''^, and the source term S to the filtering of oscillations in lower order terms in u^''^. The quadratic 
term Q is defined as follows (see [4j, Proposition 6.1 and [E]): for wi,W2 & HnH^{T), Q{wi,W2) is the 
weak limit as e ~> of the quantity 

^ exp ^"^^ ^ ^cxp ^— -L^ wi ■ Vcxp ^— -L^ W2 

+ i exp ^ fexp ( — -Lj W2 ■ Vexp (~-Lj wi 



Hence 



iV™, (1.6) 



where the resonant set Km is defined for to G Z"^ \ {0}, by 

/C„, I (fc, e Z^ J'^ + ^!^^ = ^'2' and Br, e {-1, l}^, ry^Aja + mh = msj 



and the coefficient akj.m by 

- I {{Nm, (Nk ■ V)Ni) + {Nra, {Ni ■ V)7Vfc)) . 

In order that the equation on is defined unambiguously, the value of the source term S has to be 
specified. In the present case, we have 

S = -^SB{ul)-,ypSTia), 

where Sb :'H ^Ti. is a linear continuous non-negative operator (see [4l[5l[T7]) recalled in formula (|3.10p 
below, and S't(o') depends on the time oscillations in the wind-stress a. Thus, in the next paragraph, we 
precise the assumptions on the wind-stress a'^ , and we define the source term St- In the above formula 
and throughout the article, the subscripts B and T refer to top and bottom, respectively. 

1.2. Definition of the limit equation. Let us first introduce the hypotheses on the dependance of 
the wind velocity cr^ with respect to the time variable. Since the Coriolis operator generates oscillations 
at frequencies of order , it seems natural to consider functions a'^ which depend on the fast time 
variable t/e. The case where this dependance is periodic, or almost periodic, has been investigated by N. 
Masmoudi in [T^ in the non-resonant case, that is, when the frequencies of the wind-stress are different 
from ±1. The results of [17] were then extended by the author and Laure Saint-Raymond in [5]. In 
fact, it is proved in [5] that when the wind-stress oscillates with the same frequency as the rotation of 
the Earth (i.e. ±1), the typical size of the boundary layers is much larger than the one of the classical 
Ekman layers. Moreover, a resonant forcing may overall destabilize the whole fiuid for large times. Here, 
we wish to avoid these singular behaviours, and thus to consider a more general non-resonant setting. 

Let {E,A, mo) be a probability space, and let {Or)TeR be a measure preserving group transformation 
acting on E. We assume that the function can be written 

a''{t,Xh) ^ (J (t,-,Xh:,uj] , i > 0, Xh GT'^, lo G E, 



£ 

and that the function a is stationary, i.e. 

(T(i, r 4- s, Xh; uj) = (j{t, r, Xh] OsUj) 

almost everywhere. 
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The periodic setting can be embedded in the stationary (ergodic) setting in the following way (see 
[19]): take E = R/Z ~ [0, 1), and let toq be the Lebesgue measure on E. Define the group transformation 

{Or)res. by 

ers = s + T mod Z V(r, s) eRx E. 

Then it is easily checked that 9t preserves the measure toq for all r S R. Thus the periodic setting 
is a particular case of the stationary setting; the almost periodic setting can also be embedded in the 
stationary setting, but the construction is more involved, and we refer the interested reader to [I9]- 

The interest of the stationary setting, in addition to its generalization of the almost periodic one, lies 
in the introduction of some randomness in equation (|1.2p . Hence, we also expect to recover a random 
function in the hmit e,i^ ^ 0. In fact, we wih prove rigorously a strong convergence result of this kind; 
additionally, we will characterize the average behaviour of u'^''^ in the hmit. Thus, one of the secondary 
goals of this paper is to derive some averaging techniques adapted to highly rotating fluids, which may 
be of interest in the framework of a mathematical theory of weak turbulence. 

Since the function a is not an almost periodic function, we now introduce a notion of approximate 
spectral decomposition of a. For a > 0, deflne the operator by 

Taa : A e R -!- / exp{-a\T\)e-'^^ a{T) dr, (1.7) 

and define the family of functions {aa)a>o by the formula 

CTa(T) J exp(-a|A|)e'^^J^„cr(A) dA. (1.8) 

It is proved in Appendix A (see Lemma A.l) that the family {aa)a>o converges towards cr, as a 0, 
in L°°([0,Ti] X [0,T2] x iJ, L^(T^)) for all Ti,T2 > 0. In order to simplify the presentation, we assume 
from now on that a only has a finite number of horizontal Fourier modes. This is not crucial in 
the convergence proof, but the non-resonance conditions on a are somewhat simpler to formulate in this 
case. We refer to Remark [Ol for more general assumptions. 

From now on, we assume that for all a > 0, T > 0, the function J^qCT belongs to L°°{[0,T] x E x 
T^, n L°°{R\)), and that the following non-resonance hypotheses hold: 
(HI) For all T > 0, 

VT > 0, sup ||J^qct||l^([o,t]xSxt2,li(Ra)) < +oo- 
(H2) There exist neighbourhoods V± of ±1, independent of a > 0, such that 

VT > 0, lim sup |l.?^aCr(A)||Loo([o,T]xExT2) 0. 

We refer to Remark 11.31 below for some details about the meaning of hypotheses (H1)-(H2) for 
almost periodic functions. The interested reader may also consult [13] for a treatment a resonance 
phenomena for functions with a continuous spectrum; notice however that the context in [13] is somewhat 
different, since it deals with functions whose Fourier transform is well-defined. 

Let us now explain how random oscillations are filtered: 

Proposition 1.1. Let <j> e (Rt , L'^ (E)) be stationary, and let A e R. Then the family 

(j)^ : uj e E ^ [ 0(r, uj)e-'^^dT, 0>O 
^ Jo 

converges, almost surely and in L'^{E), towards a function denoted by ^^a[0] G L'^i.E) as 9 ^ oo . More- 
over, satisfies the identity 

almost surely in uj, for all t eR. 
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Additionally, if a satisfies (H1)-(H2), then 



£x[<j] = (1.9) 

for X in a neighbourhood of ±1. 

Proposition 11.11 is proved in Appendix B, except property (|1.9p . which will be proved in the course 
of the proof page [26l 

With the above definition of £\, the source term St is defined by 

y/aaia2 ^^^^ 



where 



a{t,T,kh\uo) = I <j{t,Xh;uj)e '''''^'""^ dxh- 

J2 



Notice that S't(o') is a random function in general, and is well-defined in L^^{[0, oo) x E, L^(T)) thanks 
to (H1)-(H2) provided cr e L°°([0, T] x [0, oo) x x E) for all T > 0. 

• We now state an existence result for the limit system, based on the analysis in [4j. To that end, 
we introduce the anisotropic Sobolev spaces H^'" by 

H'''^' := {u £ L\T), \fa e N^, < s, jagl < s' , V^'^S^^u G L\T)} . 

Then the following result holds: 

Proposition 1.2. Let iy,e,l3 > be arbitrary. 

Let uo € H n and let a e L'^^{[0, cx))t, L°°([0, oo)^ x x E)). 

Assume that the hypotheses (H1)-(H2) hold. 

Then S't(o') € LJ^^([0, oo)f, iJ"'^)), and consequently, the equation 

dtw + Q{w, w) - AhW + ^J^Ssiw) + vfiSricr) = 0, ^^^^ 

W|t=0 = Mo 

has a unique solutionw e L°°(£;, C([0, oo), 7YniJ°'^)) such thatVhU belongs to L°°{E, Lf^^{[0,oo), H°''^)). 

Remark 1.1. (i) Notice that the function w is random in general because of the source term St. 
(ii) Ln J^, Proposition \1.2\ is proved for St^ — (see Proposition 6.5 p. 145). As stressed by the authors, 
the result is non trivial since the system (|1.10p is similar to a three-dimensional Navier-Stokes equation, 
with a vanishing vertical viscosity. The proof relies on two arguments: first, a careful analysis of the 
structure of the quadratic term Q shows that the limit equation is in fact close to a two-dimensional one. 
Second, the divergence-free property enables one to recover estimates on the vertical derivatives on the 
third component of the velocity field, and thus to bypass the difficulties caused by the lack of smoothing 
in the vertical direction. 

In fact, the proof of Proposition can easily be adapted from the one of Proposition 6. 5 in JJj, 
and is therefore left to the reader. The method remains exactly the same, the only difference being the 
presence of the source term St in the energy estimates. This does not rise any particular difficulty, 
thanks to the assumptions on a. 

1.3. Convergence result. 

Theorem 1.3. Assume that v = 0{e), and that ^/evfi = 0(1). 

Lei cr e iyi'°°([0,oo)^ X [0,oo)t,L°°(T2 x E)) such that (H1)-(H2) are satisfied. 

Letu^ " e L°°(i?, C([0,oo),L^)nL^^^([0, oo),/f^)) be a weak solution of l|1.2p . supplemented with the 
conditions (|1.3p and the initial data uj^'^p — uq £ Ti x H^'^. Let w be the solution of (|1.10p . Then for 
all T>0, 

u' " -ejip(-^Ljw 
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in £2([0,T] X E,H^^°)r]L°°{[0,T],L^{E x T)). 

In the case of the nonresonant torus (see (|1.12p below), it is likely that the hypothesis f = 0{e) 
can be relaxed. Indeed, in this case, the equation on w decouples between a nonlinear equation on the 
vertical average of w on the one hand, and a linear equation on the vertical modes of w on the other (see 
paragraph 11.41 below, together with Section (H). Moreover, it can be proved that the purely horizontal 
modes of w decay exponentially in time at a rate exTp{—ty/iy/e), and the rate of decay does not depend 
on the particular horizontal mode considered. Thus, in this particular case, the regime z/ ^ e may be 
investigated, using arguments similar to those developed in [5]. 

Remark 1.2. In fact, the above Theorem remains true even when the number of horizontal Fourier 
modes of a is infinite. In this case, it can be checked that the non-resonance assumptions (H1)-(H2) 
have to be replaced by the following: there exists s > 4 such that 
(HI') For alla>0,T> 0, T^a S L°°{[0,T] x E, L^{Rx, H^iT"^))), and 

VT > 0, sup \\^a<^\\L'^([O.T]xE,L^{R,H-(T^))) < +00. 

Q>0 

(H2') There exist neighbourhoods V± of±l, independent of a> 0, such that 

VT > 0, lim sup |l.FQ(T(A)||ioo([o,T]x_E._f/-(T2)) = 0. 

The H"^ regularity assumption stems from the regularity required in the stopping Lemma A. 2 in the 
Appendix. 

Furthermore, the regularity assumptions on the function a become 

fxe LJ^,([0,oo)t,L°°([0,oo), X S,i/-V2(T2)), dr<J e H\T^,L°°{[0,oo)t x [0,oo), x E)). 

Remark 1.3. Let us now explain the meaning of hypotheses (H1)-(H2) for almost periodic func- 
tions. Let kh e I?, and let G L°°([0,cx)) x T^) such that 

</.(T,x„) = e<--'' ^ </>^e'''^ 

where M is a countable set. The fact that 4> as only one horizontal Fourier mode is not crucial, but 
merely helps focusing on the time spectrum. Then it can be checked easily that for all a > 0, 

2^ a^ + it^- A)2 

In particular, there exists a constant C > such that 

||-?^a0||L°°(T^Li(R;^)) < ^ \(l)f,\ 

Thus hypothesis (HI) is satisfied provided X^/^eA/ I'^pl < 
On the other hand, assume that 

?7:=d(M,{-l,l})>0, (1.11) 

i.e. that there are no frequencies in a neighbourhood of ±1. Then if \ <E (—1 — ri/2, —1 + 77/2) U (1 — 
77/2, 1 + 77/2), we have 

|A-m|>| V/iGM, 
and consequently, setting :— (±1 — 7]/2. ±1 + t//2), we have, 

sup \\Ta4>{\)\\Locn2-.<C-. 

\ev-uv+ V 
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Thus hypothesis (|l.lip entails (H2). Additionally, hypothesis cannot be easily relaxed, as shows 

the following construction: consider the sequence ^„ :— 1 — 1/n, and choose a sequence of positive numbers 
(bn such that 



< oo. 



For r e M, set 



Then for all a > 0, for all k > 



ir) 



^a?!'(/i/c) = — TV — 7T2 ^ 



a 



In particular, 

lim J^aHf^k) — +00 

for all k, and thus condition (H2) is not satisfied. 

1.4. Average behaviour at the limit. We have already stressed that the solution w of equation 
(|1.10p is, in general, a random function. Thus one may wonder whether the average behaviour of w at 
the limit can be characterized. In general, the nonlinearity of equation l|1.10p prevents us from deriving 
an equation, or a system of equations, on the expectation of w, which we denote by E[w]. However, when 
the torus is non resonant, equation (|1.10p decouples, and in this case we are able to exhibit a system of 
equations satisfied by E[w]. 

Let us first recall a few definitions: 

Definition 1.4 (Non-resonant torus). The torus :— x [—a, a) is said to be non-resonant if 
the following property holds: for all (k, n) £'L^\ {0} x \ {0}, 

(377 e {-l,l}^7?lAfe +?72A„-fe - ?73A„ = O) ^ fcaJ^sl^s - h) = 0. (1.12) 



We refer to [3j for a discussion of hypothesis l|1.12p and its consequences. Let us mention that (|1.12p 
holds for almost all values of (ai,a2,a) € (0,oo)^. When the torus is non-resonant, the structure of the 
quadratic form Q defined by l|1.6p is particularly simple, and the system (|1.10p can be decoupled into 
a two-dimensional Navier-Stokes equation on the vertical average of w, and a linear equation on the 
z-dependent part (see [4]). The advantage of this decomposition in our case is that the vertical average 
of St{<7) is deterministic, at least when the group transformation (0r)T>o acting on E is ergodic (see 

Definition 1.5 (Ergodic transformation group). Let (6'i-)7.gR be a group of invariant transforma- 
tions acting on the probability space {E,A,mo). The group is said to be ergodic if for all A £ A, 

{0rAcA Vr e M) ^ TOo(A) = or TOo(A) = 1. 

We now state the result on the average behaviour at the limit: 
Proposition 1.6. Assume that the transformation group (0t)tgb is ergodic. 
LetuQ e TinH^'^, and let a € L°°{[0,oo)t, x[0,oo)r x E x T^) such that the hypotheses of Theorem fTJI 
are satisfied. Let w € L°°{E, C{[0, oo), Ti. H'^'^) Lf^^{[0, co), H^''^)) be the unique solution of equation 

(fTTnli . 

Let w = (w/^,0) e C([0,oo),L2(T2)) n L]^^[[Q,oo),H^{T^)) be the solution of the 2D-Navier-Stokes 
equation 

dtWh + Wh ■ VhWh - AhWh H -j=J -Wh + vIM [5't(ct)]^ = V^p, 

av2 V £ 

AivhWh = 0, 

1 r 

Wh\t=o{xh) - / ua,h{xh,z) dz. 

a ./n 
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Then the following properties hold: 

1. As e,v ^ as in Theorem ] 1.3[ we have 



w in L^{[Q,T]xT X E). 



In particular, the weak limit of u^'^ is a deterministic 
2. Assume additionally that the torus T'^ is non resonant. Then 

V\w\ = w + w, 

where w solves a linear deterministic equation 

dtw + 2Q{w,w) ~ AhW + ^Sb{w) = 0, 
w\t=o = Uo~ W\t=o- 

1.5. Strategy of proof of Theorem II. 3L The proof relies on the construction of an approximate 
solution, obtained as the sum of some interior terms - the largest of which is eiq){—T L / e)w{t) - and some 
boundary layer terms which restore the horizontal boundary conditions violated by the interior terms. 
We refer to the works by N. Masmoudi [El [17], N. Masmoudi and E. Grenier pj], N. Masmoudi and 
F. Rousset [E], and F. Rousset [22j for an extensive study of boundary layers in rotating fluids, or in 
incompressible fluids with vanishing vertical viscosity for . We emphasize that in fine, all terms will 
be small in norm, except exp{—TL/e)w{t). 

Following [4] (Chapter 7), let us assume that as e, ^ 0, 

p-.;^ipint + ipBL^pi„t,0^ (1-13) 



e £ 



where 



BT/ f ^ a — z\ ft z\ 

U (t, Xh,Z) = UT [t,-,Xh, ] + UB [t,-,Xh, - ] , 

BT/ ft a — z\ ft z 

p [t,Xh,z) ^pt t,-,Xh, \+Pb\ t,-,Xh,- 

\ £ JJ J \ £ 77 

Above, 77 is a small parameter that will be chosen later on. The function UT{t,T,XhX) is assumed to 
vanish as C ^ 00 (same for pt,Pb, ub)- 

We then plug the Ansatz (|1.13p into equation l|1.2p . and identify the different powers of £. In general, 
there is a coupling between u'"* and u^^: indeed, we have seen that it is natural to expect that 

U{t,T) = e^p{-TL)w{t), 

at first order, and thus m'"* does not match the horizontal boundary conditions in general. As a conse- 
quence, the value of m'"* at the boundary has to be taken into account when constructing the boundary 
layer term vP^. On the other hand, because of the divergence-free constraint, the third component of 
u^^ does not vanish at the boundary, which means that a small amount of fluid may enter or leave the 
interior of the domain. This phenomenon is called Ekman suction, and gives rise to a source term 
(called the Ekman pumping term) in the equation satisfied by w'"*. This leads to some sort of "loop" 
construction, in which the boundary layer and interior terms are constructed one after the other. 

The first step of this construction lies in the definition of boundary layer terms. In the periodic case, 
this is well-understood (see |4l[l7l[5]); thus the main contribution of this article in this regard lies in the 
definition of boundary layers in the random stationary case. Hence, Section [2] is entirely devoted to that 
topic. Section[3]is concerned with the definition of first and second order interior terms; in particular, we 
derive in paragraph 13.21 the limit equation for the system (|1.2p . In Section [H we prove the convergence 
result, after defining some additional corrector terms. At last, we prove Proposition II .61 in Sectional 
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2. Construction of random boundary layer terms. The goal of this section is to construct 
approximate solutions of equation (|1.2|) . which satisfy the horizontal boundary condition at z = a, and 
which are localized in the vicinity of the surface. Such a construction has already been achieved in the 
case when the function a is quasi-periodic with respect to the fast time variable (see [iTj in the non- 
resonant case, and [5] in the resonant case). Thus our goal is to extend this construction to a random 
forcing. The main result of this section is the following: 

Lemma 2.1. Assume that P^/eV ^ 0{\). Le^ cr e M^i'°°([0, oo)t x [0, cx))^, L°°(T2 x £;)) he such that 
(HI )-(H2) are satisfied. Then for all 5 > 0, there exists a function u!^^'^ £ L°°{[0, oo)t xT x E) which 
satisfies the system 

OtUrj. + -63 A Uj, ' - vd^Uj, - AhUj, =CII1 + -I [eiy) ' \\a\\ I 



L°°([0,oo)x£;,L2(x)) 



BL.S o 

Ayvu^^'^ = 0, 



and such that 

II BL.Sw ^ I \ — 1/411 -Bi,i5|i _ I 

supljuy ■ II 00 < 00, sup(ei^) ' \\uj, ||l~([o,oo)x£;,l2(t)) < +00. 

<5>0 <5>0 

Moreover, u^^^^^ is exponentially small. 

The above Lemma entails in particular that for all S > 0, u^^'* is an approximate solution of (|1.2p . 
which satisfies the appropriate horizontal boundary condition at z = a. The Lemma is proved in the two 
next paragraphs: we first explain how u^^'^ is defined, and then we derive the and L°° estimates. 

2.1. Construction of the boundary term at the surface. As explained in the Introduction, 
the idea is to consider an Ansatz of the form 



BL,5 



. t a ~ z 

{t,Xh,Z,UJ) ^ UT [t, ~^^h, ^ ;UJ 



where is a small parameter (whose size has to be determined) and 

lim UT(i, r, a;/i,C;a;) = M t,T,Xh,Ljj. 

C, — ^cjo 

Hence we expect to be of order 77/?||ct||oo in L°° . Moreover, the divergence-free condition entails that 
the third component of ut is given by 

/•oo 

UT,3 (C) = -V <i^^hUTJiiC)dC] 



thus ut,3 = ©(/Jry^llcrllvFi.^")- At last, in order to be consistent with l|1.13p . we assume that the pressure 
inside the boundary layer is given by 

1 

p[t,Xh,z,uj) fa -pT 



where pr = 0(/3ry| |cr| |oo)- Then the pressure term in the third component of (|1.2p is of order /3||cr||oo/£, 
whereas the lowest order term in the left-hand side is of order ?7^/3| |cr| |vyi,oo /s. Thus, since ry is small, we 
infer 

d(;PT = 0, 

and since pt vanishes at infinity, we have pt = 0: at first order, the pressure does not vary in the 
boundary layer. Thus, we now focus on the horizontal component of ut, which is a solution of 

dr - + f-""'^) = 0, (2.1) 

dcUT,h\c=o = -vP<^{t, X, y, uj), (2.2) 
utji\c=+oc. = 0. (2.3) 



IQ 



We now choose rj so that all the terms in l|2.ip are of the same order, that is, 

77 = 

Moreover, since cr is a stationary function of time, it seems natural to look for stationary solutions of 
(|2.ip . and thus for fundamental solutions ipi, ip2 of (|2.ip in the following sense: ipi (i — 1, 2) is a solution 
of Il2.ll) in the sense of distributions and satisfies 112.31) , and 





where 5q denotes the Dirac mass at r = 0. If we can construct ipi and (^2 satisfying the above conditions, 
then a good candidate for ut is 

noo 

UT,h{t,T,XhX\^) ^ -VveP X! / <yj{t,T - s,Xh;uj)(pj{s,C)ds. 

ie{i.2}-^o 

Hence we now define ipi,ip2. Since the fundamental solution of the heat equation is known, let us make 
the following change of unknow function (see [17]): 



J = 1,2. 



Then, setting ef := (l,=Fz), := (±i, 1), we infer that — Ge^ , where G satisfies 

drG-d'^G^O, T>0,C>0, 
G|c=o(t) - <5o(t), 
G|(^^_i_oo = 0. 

The boundary condition at C = should be understood as follows: for all (p E Cb(K.), for all r > 



(2.4) 



lim 



(P{t — s)G{s, Qds 



ip{T). 



It can be checked (see Chapter 4, section 1 in [T2j) that 



G(t,C) := 



C 



47rr3/2 



exp ( -|- ) for r > 0, C > 0, 



is a solution of (12.41). which leads to 



G(r,C) \e-"e+ + e^"e^ 



Unfortunately, when we integrate this formula with respect to C in order to obtain an explicit expression 
for Ut h , the convolution kernel thus obtained is 



<Pj(t, C) 



1 



/47rT 



exp 



At 



and is not integrable near r — +00. Hence, in the spirit of [TT], we consider an approximate corrector in 
the boundary layer: for (5 > 0, we set 



G5(t,C) = 



C 



47rr3/2 
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Then the corresponding corrector is given by 

roc 

Ut hit^T^Xh,C,^) ^ -PVe^ y2 / ipjis,C)exp{-6s)aj{t,T - s,Xh;uj) ds (2.5) 



je{i,2}- 



The approximate corrector satisfies the exact boundary conditions at C = 0, and equation p.ip up 
to an error term of order S 

The third component of is then given by 

/oo 
div,,4,/i> 

which yields 

where Lp is defined by </5'(C) = exp (~ ' = 0- 

In horizontal Fourier variables, we have 

43ft,r,x;.,C^)= E E^'"^"" r^f4V*(^,r-s,fc,.,c.)e-^^±-ds (2.6) 

V47raia2 ^,^2 ^ Jo Vv«/ 

where 

<7±(fc,o = *^;-^(fch)±(fc;)^-<7(fc,o- 

Now, set 

BL,5/, X s f, ^ a — z 

Urp (t,Xh, z;uj) := Uj. [ t, -,Xh, ^ ; 



e 



It can be readily checked that 



\{dt+'-dr-A,-ld'^+s)u^^^,{t,l,x,,^; 



a BL,S , 1 . BL,S a2 BL.5 a BL,5 

OtUrp H 63 AUrp — lya^Urp ■ — lAhUrp 



There remains to evaluate in L°° and L^. 

2.2. Continuity estimates. This paragraph is devoted to the proof of the following Proposition: 
Proposition 2.2. Assume that a e L°^{E x [0,oo) x T'^,Cb{^r)), and that a satisfies (H1)-(H2). 
Then for all T > 0, there exists a constant Ct > 0, such that for all S, z/, e, /3 > 0, 

114, 9c4, C9c4|Loo([o,T]xR.xT=x[0,oo),xiJ) ^ ^tV^P, (2.7) 
T' ^CU^T, C9cWt|Loo([o,t]xR.x£;,L2([o^oo),xT2)) ^ <^TVe'^/3. (2.8) 



Remark 2.1. With the assumptions of Theorem \1.3\. the same bounds also hold for all the derivatives 
of ut,s with respect to the macroscopic time variable t and the horizontal space variable x^. 

Proof. We focus on the horizontal component of u^; the vertical one is treated with similar arguments. 
Recall that ^ is given by l|2.5p ; in order to simplify we set ± ia . 
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First, we write 

' ^?r;^-(-£-^')"*<^''--""'* 



(2.9) 



(2.10) 



The term l|2.10p can easily be evaluated thanks to Lemma A.l in the Appendix; notice that since the 
convergence given in Lemma A.l is not uniform with respect to r e [0, cx^), we cannot derive an estimate 
in L°°([0, oo)t-) at this stage. Hence we keep the variable t for the time being; there exists a constant 
C > such that for all r > 0, i? > 0, 



exp 



4s 



, r - s, •, w)e'''*^ ds 



L°°([0,T]x£;,L°°(T2)) 



< C\\a — (Ja\\L°°{[0,T]xEx[T-Fi.,T]xT^) 

1 



—= cxp (5s ) ds 



4s 



(2.11) 



C\\a\ 
C,, 



L°°([0,T]xKtXT2xB) 



< yll^ ~ fa||L°°([0,T]x£;x[T-fl,r]xT2) 



+ ^ II'''IIl°°([0,T]xK^xT2x£;) 

Choosing R = 5^"^ , we deduce that 



]xT: 

exp(— (5i?) 



—j= exp \ <5s ) ds 



^exp ( -(^s ) (ct± (•,t- s,-,cj)e±'^ ds 



L°°([0,T]x£;,L=°(T2)) 



- '^"^^ '^"llL°°([0,T]xBx[T-^,r]xT2) 



c 



exp I 



As for the term l|2.9p , recalling the definition of a a , we have 



^fs 



As 



—= exp — 5s \ TaU {-^ T — s, ■,'^)e ds 



(2.12) 



We first evaluate 



^ J -i=e-"l^lexp(^-£-(5s^ J^„a±(-,A,-,a;)e*^(^-")e±^'^dAds. (2.13) 

We split the integral into two parts, one going from s = to s = 1, and the other from s = 1 to s — oo. 
It is obvious that for all C > 0, (5 > 0, A € R, 



1 1 



: exp 



4s 



< 



1 



/i'^^Py 4s 



^'^ds<l 
- 2 



(2.14) 



Integrating by parts the second integral, we obtain 



1 



As 



^expf-^)e-(^+^(^±i)> ds 



1 



6 + i{X±l) 
1 



exp I -- 



2(<5 + z(A±l))7i s 



2s 



exp 



4 s 



(2.15) 
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We are now ready to derive the L°° estimate; the function 



is bounded on M. Hence, gathering (|2.14p and l|2.15p . we deduce that there exists a constant C such that 
for all C > 0, (5 > 0, A e M, 



< C 



n — eXD f~ — ^ g-(5+i(A±l))s J 

Jo 4s 
Inserting this inequality in i|2.13|) . we obtain 

/ —j= exp [ — (5s I Ta(T^{-,T — s, •, ix')e^*'* ds 

Jo Vs \ 4s J 



1 



< C I e-"l^l 

< C 



1 



1 



\S + i{XTl)\ 

sup ||(T+_Q||ioo(£;xT2,Ll(R)) - 



J?aO-='=(-, A, -,0;) 



L«=(T2) 



|5 + i(A±l)| 



L°°(T2) 

dA 



1-5 



^-^^^||^„a(,A,,a;)||^^(,.) 



dA 



+ C / \\Ta(7{-,X, -,1^)11 ^^,j.2) dX 

Jr\v± 



< C 



sup ||J^Q(T||ioo(£;xT2.Li(K)) + sup ||J^QCr(A)||^^(T2) 



Above, we have used the following facts: there exists a constant Ci > such that 

|(5 + i(AT 1)1 > |At 1| > ci VAeM\Kt, 
and there exists another constant C2 > such that 



1 



< 



±H-C2 



1 



1(5 + i(A T 1)1 - y±i-c. v/^M^TTAy 
We deduce that for all (3; > 0, for all (5 > 0, r > 0, 



zdA < Cln((5). 



lL~([0,T]xT2x[0,oo)Cx£;) 



. I exp(-|) 
<5 



|cr - (TQ||ioc([o_T]x[r-5-i,r]x£;xT2) + sup ||J^q(t(A)|1 ln((5) 



Taking the infimum with respect to a of the right-hand side, with (5 > fixed, we deduce that 



SUp||MT,?i|Loo([o,T]x[0,oo)rxT2x[0,co)Cx£;) - C'VEJ^/?- 

We now turn to the derivation of the I? estimate, which is similar to the above computations. The 
main difference lies in the fact that we need to integrate by parts l|2.15p yet another time, which yields 



4 exp 



4s 



1 



5 + i(A± 1) 



S'^P — T 



2(5 + i(A± 1))2 



exp I -- 



2((5 + i 



14 



where 



Consequently, remembering (|2.14p . we have 



1 



+ 



Is " \ 4s 
1 

2|<5 + i(A± 1)|2 
1 



2|J + i(A±l)|^ ji S2 



1 



exp 1 



ds. 



Plugging this estimate into (|2.13p and using (H1)-(H2), we infer that for all C > 0, for all s > 0, 



1 



< c 



— exp \--^ <5s ) (Tq (•, r - s, •, u;)e ds 



L°°(T2) 



+Cexp 



4 / l^^^AGK ll-^"'^('^)lli'°([o,T]x£;,L-(T2)))ln('5) 



2 



1 

5" 

1 52 



ds 



1 + ^^^^ ll-^a'^('^)llL^([0,T]x_B,L°=(T2) 



1 



= ([0.T]x£;,L~(T2))) J I • 



Here, we have used the inequality 



d\ 



< 



T1 + C2 



d\ 



c 

< — . 



There only remains to prove that each term of the right-hand side has a finite norm. First, thanks to 
Jensen's inequality, we have 



f°° _rnl 

< 2 / ds e ^ dx < oo. 
Jo Jo 



Similarly, 



C30 / /'OC 



1 

1 S2 



OO /"OO 



ds] dC<C 



Ji S2 



dsd(^ 



\(j>{x)\^ dx\ < OO 
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We also have to evaluate the norm of the integral in l|2.1ip ; we have 





/•oo 




Jo 


< 




=Vsc, 


¥j 


u—Ss 




< 


J 


< 


5i J 


< 


c 
7a ■ 



1 / 

—= exp — Ss ] ds 



Vs 



As 



: exp 



4ii 



u du 



dx 



2 Jo Jo 



1 fx 
— exp I — u 



2u 



du dx 



X 



: exp ( — - — u ] du dx 



Gathering all the terms, we obtain, for all a,S > 0, for all r > 0, 



Khir)\ 



L~([0,T]x£;,L2([o,oo)<;,L~(T2))) 



, If ~ O"a||L°°([0,T]x[r-<5-i,Tlxg,L°°(T2)) 

Op €iy 3 



""""^K + sup |^„a+(A)| (l+\n{S) 
d2 \eV- \o 



Taking the infimum of the above inequality with respect to a, we infer the estimate on u^/j. The 
estimates on u^^, d^u^, and Qd^u^ are derived in a similar fashion. □ 

Remark 2.2. Stationary boundary layer terms relative to Dirichlet boundary conditions can also be 
defined: consider for instance the boundary condition 



The construction is the same as for Neumann boundary conditions, and is in fact more simple because vie 
need not integrate with respect to the variable (. Thus, with the same notations as above, the boundary 
layer term at the bottom is given by 



1 

h{'t,T,Xh,C,<^) = 2 X! / Gs{s,Q[e~'''ej +e+'''ej]cB,j{t,T- s,x,y;uj)ds, 



ie{i.2} 



and 



ve 



EE 







exp 



4s 



-5s±is 



ds. 



2.3. Previous results in the quasi-periodic case. For the reader's convenience, we have gath- 
ered here previous results appearing in ^ [l^ , in which the authors compute the boundary layer term at 
the bottom of the fluid. We recall that it can be expected that the solution in the interior of the domain 
behaves like some function exp(— with w € L°°(£',C([0, oo), 7i)). In general, the horizontal 
component of such a function does not vanish at z = 0, and thus a boundary layer has to be added in 
order to restore the Dirichlet boundary condition. Consequently, we seek for a boundary layer term 
which is an approximate solution of equation l|1.2p and which satisfies 

<HU=oii^'"^)= E CB,h{t^k)e<--'^e-^>^>'"'; (2.16) 
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for the boundary layer term at the first order, the coefficient CB.h{t, k) will be given by the formula 

However, we will also use this construction for the lower order boundary layer terms, and thus we keep 
an arbitrary value for CB,h{t, k) for the time being. 
As before, we assume that 

Xh, z) = UB (t, -,Xh, 



The decomposition (|2.16p leads us to search for a corrector ub satisfying 
where each term UB,h,k is an approximate solution of (|1.2p . and 

The periodicity in time of the boundary condition prompts us to choose UB,h,k as a periodic function of 
r, with frequency Afc. Also, it is classical to seek UB,h,k as an exponentially decaying function of C,; the 
rate of decay is then dictated by the equation. The precise expression oiuB,h,k is the following (see [T7]): 
• First case: kh 0. 

In this case, UB,h,k is an exact solution of (|2.ip . and is equal to 

UBMAt, r, X, y,0=Yl ^fc cj)e-'^^^+^<-^^-^i'^ (2.17) 



vt = VlT Afc 



l±i 

"71"' 



^±j.^.^-j _ 1 fcBAt^k) ±icBMi^k)\ _ CB,i{t,k)±icB,2{*^ k) ( 1 



where 



2\cBAt^k)^icB,i{t,k)) 2 \Ti 

The vertical part of the boundary layer is then given by 

UB,z,k{t, T, X, = Vi^E ^< ■ < u)e~'^>^-^+<----^'i^. (2.18) 

± Vk 

• Second case: kh = 0. 

In this case, the construction of the resonant boundary layers in ^ proves that there are indeed boundary 
layers, but which are of order and not y/ei' in general. Thus the size of the boundary layer depends 
(slowly) on time. 

First, notice that if kh — 0, then — — sgn(fc3) — ±1. As in the first case, we decompose CB.h{t, k) 
onto the basis (l,±i) : 



cnAt^ ^) = \Y1 (^B,i{t, k) =F icB,2{t, k)) . 

As a consequence, we have 

E CBAt,^.k^)e''^-^ = a+{t)e'^ (j) +a_(t)e-^^ 
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where 



a±it) = ^ (cB,i(t, 0,^3) =F»CB, 2(^,0,^:3)) 

fc3,sgn(fc3) = ±l 

7±W= {cB,l{t,0,k3)±lCB,2it,0,k3))^ 



fe3,sgn(fc3) = ±l 



The terms 7±e (1, ^i) give rise to a classical boundary layer term, namely 



E7±We±— "^^(^J, with,7± = l±z. 



± 



For the terms Q;±e^''^(l, ±i), we rather use the following Ansatz (see fS]) 

u^^''^^%t,x,,,z) = ^ f')j2a±{t)e^'^ (}) . (2.19) 



In order that u^^''"^^ is an approximate solution of (the linear part of) equation (|1.2p . the function V' 
must be such that 

-|^'(X)-^"(X) = 0, 

which yields 



With this definition, u^^''"^^{t) vanishes outside a layer of size \/1a localized near the bottom of the fiuid. 

B 



Hence M^L.res jg approximate solution of the linear part of equation (|1.2p . and is exponentially 



small. 

Now, set 



fcez^ ± ^ ■ 



The complete boundary layer term at the bottom is given by 



U^^{t,Xh,z) =UB ( t,-,Xh,^= ] +U^^'''''\t,Xh,z) 



We now give some estimates on the boundary layer terms constructed in this paragraph: 

Lemma 2.3. Let ub he defined by (|2.17p - (|2.18p and u***** by l|2.19p . Then the following estimates 
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hold, for all t > 



( 



\uB,h(t)^ C9cMS,/i(i)||ioc([o,oo),,L2(T2x[0,oo)<;)) 



< c 



\feh#0 



\uB,h{t), C9c"i3J>WllL-([0,oo).xT2x[0,oo),) ^ ^ \cB^h{t,k)\ 



( 



5] |cB.,(t,fc)|2^|fc3 



and 



,BL,res 



The proof of the above Lemma is left to the reader. Notice that according to the definition of rj^, 
we have 

C^<|r/±|<1 VfceZ^. 
|fc| I I 

Corollary 2.4. Assume that there exists N > such that 

CB,h{t,k)^Q if\k\>N, yt>0. 

Then the boundary layer term is an approximate solution of the linear part of equation (|1.2p . Pre- 
cisely, there exists a constant Cn , depending only on N , such that 



L°°([0,T]X2(x)) 

<Cnv^I^ sup ( V (|cB.„(t,fc)|2 + |a,cs./.(t,A:)|2) 



1/2 



Proof. By construction, is an approximate solution of the linear part of equation l|1.2p . with an 
error term equal to 



[{dt - A,i)ub] ( t, Kxh, -j^ ) + <^ 



where dt is the derivation operator with respect to the macroscopic time variable. Thanks to the as- 
sumption on the coefficients cbji, we have 



/o JtJe 



dtUB { t, - ,Xh, 



e Jev 



-I 1/2 



dmo^Lu) dz dxfi dt 



<Cn 
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whereas the term AhUB is bounded in L°° {E, L'^{[0,T], H 
term due to u^L.res satisfies 



')) by CNiei^y/^CBMl At last, the error 



hCBM 



L°=(£;,L2(T)) 



3. The solution in the interior at main order. This section is devoted to the construction 
of the first order interior terms in expansion (|1.13p . At this stage, we merely know how to construct 
boundary layer terms which deal with the horizontal part of the boundary conditions l|1.3p : moreover, 
following the analysis in paragraph [TTTl we expect u'^'^{t) to behave like eyi-g{~tL / £)w{t) in the interior. 
Thus the idea is to define a function 



:=exp 



t 

t,- 

e 



where w'"* and Su™^ are corrector terms, such that u'"* is an approximate solution of (|1.2p . We also 
require that 



BL,(5 



,BL 



satisfies the boundary conditions l|1.3p at main order. 

Let us now explain the role of the correctors v'"* and 6u™^ : it can be checked in the formulas of the 
previous section that the third components of the boundary layer terms do not vanish at the boundary: 
indeed, one has 



u. 



BL,(5 
T,3|z=a 



0(/3eHkll 



BL 



^,^OiV^\\w\\). 



The role of the corrector w'"* is precisely to lift these boundary conditions and to restore the zero- 
fiux conditions at the bottom and at the surface. Consequently, the term d'"* has fast oscillations (at 
frequencies of order e~^), and in general, w'"* is not an approximate solution of l|1.2p . Filtering out the 
oscillations in the term 



dv" 



1 



at e 



yields the source terms Sb and St in the equation satisfied by w (see equation l|5.ip ). The remaining 
oscillating terms in the expression above are then taken care of through the addition of the corrector 

The organization of this section is as follows: first, we deal with the corrector d'"*, by giving a precise 
definition and explaining how oscillations are filtered. Then we derive the equation on the function w\ 
in general, this equation depends on the small parameter (5, introduced when constructing the boundary 
layer terms at the surface. Thus in the third paragraph, we identify the limit as (5 ^ of the function 
w, which yields the envelope equation. Eventually, the fourth and last paragraph is concerned with the 
definition of the corrector 5u™^. 

3.1. Lifting the vertical boundary conditions. In the rest of this section, we set 

CB,:i{t,T,Xh) := -—j=UB,'i\C^=o{t,T,Xh) 



ik'h ■ w 



k_ ik'^-Xh p-iy^kT 



where 



wt{t)^--{NkMt)) 



ni{k) ± m2(fc) 
n2{k) q= mi(fc) 



2Q 



and 77^ was defined in the previous section. In order to shorten the notation, we set 

Ev^ ifcl • wt(t) 
2^ ^ ^J"^ ' e 



«3 

The function w will be defined in the next section. 
Similarly, we set 

1 



fcaGZ ± '^fc 



1 f°° 

= o5I / [div/iO-Tirot/iCj] (t,r-s,a::,,;w)e"''^=^'''ds. 
^ _l_ Jo 

We also set CT,3{t,T, kh;uj) = CT,3{t,Xh',uj)e~'''^'^'^'^dxh, so that 

j. Jo 

With the above definitions, the function cb,3 is quasi-periodic with respect to the fast time variable 
r, whereas ct,3 is random and stationary with respect to the fast time variable r. 

• Defintion of t;'"*. We now define a function w'"* which is divergence free and such that 



(3.1) 



Of course, conditions (|3.ip do not determine w'"* unequivocally. A possible choice is 

W3"*(t, T, x) = ^ [ei^/3cT,3 (i, T, Xh) z + y/evcB,3 {ti T, Xh) [a - z)] , (3.2) 
<*(t, T, x) = \^h^l^ W^CB,3 {t, T, Xh) - eu(3cT,3 (i, T, Xh)] . (3.3) 

In fact, since cb,3 is an almost periodic function, a more convenient choice can be made, which is the 
so-called "non- resonant" choice in [17]. In this case, the equation on 5u™^ is slightly more simple, since 
there is no source term due to cb,3- However, we have chosen here not to distinguish between stationary 
and almost periodic boundary conditions, and thus to work with the expressions l|3.2p . p.3p . 

• Filtering the oscillations. We give here the statement and proof of a Lemma which will be 
useful in the construction of 5u™^ and w. 

Lemma 3.1. Let T > be arbitrary. 

Letv£ L°°([0,r] X [0,oo)^,L2(T2 x E)) such that drV e L°°{[0,T] x [0, oo)r, L^{T^ x E)) and 

divv = 0, (3.4) 
vslz^oit, T, Xh) = \/evcB,3{t, T, Xh), (3.5) 

V3\z=a{t, T, Xh) = PeiyCT,3{t, T, Xh). (3.6) 

Then as 9 oo, the family 

Se:= \ ! C{-t)^ \drv + eg A v] dr 
9 Jo 

converges almost everywhere and in L°°([0, T], L^(T^ x [0, a] x E)), and its limit does not depend on the 
choice of the function v. Precisely, 

lim Se =: S[cb,3, ct.s] = , ^ Y ^^-a, [^CBAkh) - {-l)''^ f3evcT,3{kh)\ Nk. (3.7) 
e^oo y/aaia2 \k'Y 
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Remark 3.1. In the above Lemma, the operator £\, which was originally defined for random sta- 
tionary functions, has been extended to almost periodic functions: if 

c{t) = "A^^'"'"' 

with la^l < oo, then 

I re _ 

E\[c]:^ Van - I e ^^c(t) = lA=^a^. 

Proof Let w\ solutions of and let ^ = - w^. Notice that V € L°°{[0,T] x 

[0,oo)r; L'^{E,n)), and drV G L^{[Q,T] x [0, c»)^; x T)). We write 

C{-tW [drV + 63 A y] = C{-t) [drV + LV] ^ [C{-t)V{t)] . 

Consequently, 

^ a-r) [drV + P (ea A V)] dr = ^^''^^7'^-" . 

The right-hand side of the above equality vanishes in L°°([0, T] x E,L'^{T)) as 00. Hence the limit 
is independent of the choice of v. 

In order to complete the proof of the lemma, it is thus sufHcient to show that the limit exists for the 
choice (|3.2p - (|3.3p . and to compute the limit in this case. First, we recall that for any function F e L^(T), 
we have 

FF= J2 iNk,F)Nk. 



It can be readily checked that if kh = 0, then {Nk, 9rW'"*) = 0. Thus for all k = {kh, ks) e such that 
kh 7^ 0, we compute 

(7Vfc,9rw'"*) = - / cos{k'^z)nh{k) ■ — ^ {evl3drCT.?.{- ,kh) - \fevdrCB,i{- ,kh)) dz 
1 



+ - / "3(fc) sin(fc3z) [^/evdrCB,?,{-,kh){a - z) + ev(3drCT^?,{-,kh)z) dz 



+ lfe3=o rihik) ■ — ^ (eiypdrCT,3i-,kh) - VevdrCB,3{- ,kh)) 
\'^h\ 



Notice that if /cs = 0, then 



nh{k) ■ k'^^0; 

consequently, we have 



In a similar way. 



(Nk, 63 A 1)'"*) = - / cos(fc3z)nh(fc) • -jvHr {£vl3cT,3{--, kh) - ^/evcB.zi-, kh)) dz 
a Jo \Kf^\ 



lk3=onhik) ■ -rrrw- {'^'^PcT^ai', kh) - Ve^CB^si', kh)) 



Wh 



U3=0 1 



— {^fevcB, ?,{■■, kh) - ev(3cT,3{--, kh)) 
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We deduce from the above calculations that 



(3.8) 



H , lk^=oj^ (Vei^cB,3{t,T,kh;uj) - eiyPcT,3{t,T,kh;uj)) Nk- 

We decompose the sum in the right-hand side into two sums, one bearing on kh such that |fc;i| > A, 
denoted by Si^a, and the other on \kh\ < A, denoted by 5*2, a, for some A > arbitrary. Using the fact 
that iiyjev = C(l), we have 



l|5l.A(i,T)||i. 



< Cev 



Cev 



^ k'\k' 



\kh\>Ak'j,& 



1,= 



TTT {(^V£^CT,3{t,T, kh]uj) - CB,3{t,T,kh;Lj)) Nk^,o 



\kh\>A 



|fc„|>A 

+Cei^ ^ (|c_B,3(<, r, fc/i;tj)P + |cT,3(t, r, fc/i;w)p) . 

Since cb, ct, drCB, drCr belong to L^(T^, L°°{[0, oo) x [0, T] x i?)), we deduce that the sum Si^a vanishes 
in L°°([0, T] X [0, (X)), L^(T2 x [0, a] x E)) as A^oo. Thus we work with A sufficiently large, but fixed, 
so that Si^A is arbitrarily small in norm, and we focus on 5*2, a- 
For k fixed, we have, according to Proposition ll.il 

e'^''''drCT,3it,T,kh;uj) dr = ^^^''q J e"^"^ CT,3{t,T,kh;uj) dr 

+ \ {e'^'''cT,3{t, 0, kh\U0)- CT,3{t, 0, kh-, io)} 
o 

— * -iXk£-\k [cT,3it,kh)] (uj) 



in L°° {[0, oo) t, L'^{E)). The calculation for cb,3 is identical. 
Using Lebesgue's Theorem, we deduce that as — > oo 



i S2.A{t. r) dr — ^ V V ^f-A, [V^CB,3{t, kn) - (-l)'=^ez./3cT,3(i, A:,)] TV^, (3.9) 
and the convergence holds in L°°([0, T], L^(T x E)). Moreover, we have 

E Sl^-AjcT^3(i,fc.)]r <C Y:pj^llf-A.[cT.3(i)]|li.(,.) 



<c\\ 



A similar estimate holds for cb,3- Thus the right-hand side of (|3.9p converges in L^(T x i?) as A ^ cxd. 
Eventually, we infer (|3?7i) . □ 



• Computation of the source terms. For the sake of completeness, we now derive an expression 
of S[cb.3, ct,3] in terms of w and a. We begin with f-A^ [cb,3]. Remembering the definition of CB,3(t, ^/t), 
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we have 



Easy calculations lead to 



± '7? 



Thus, we define the Ekman pumping term at the bottom of the fluid by 

Seiw) {Nk,w)AkNk, (3.10) 

fcez^ 

where 

2^/2a|fc'|2 yr^A^ ^ ' 

There remains to compute the coefficients £^\^{cT.3{t,kh)); since the boundary condition ct.s de- 
pends on the small parameter S, the corresponding Ekman pumping term will depend on 5 as well. 
The limit as 6 vanishes of the corresponding source term will be computed in the next paragraph. By 
definition of 8\, we have, for all kh € Z^, for all A G M, 

- \ E lii- r f khMe-^^ dr^ ds. 

where 

a^{kn) = (ik'h±{k'h)^)-a{kn). 
Thanks to Lebesgue's dominated convergence Theorem and Proposition I we infer, for all (5 > 0, 

1 r°° 

£x [cT,3(t,fc/.)] (^) - o E / [^^(.t^kh)] ie^si^)e-''^'' ds 
^ _i_ Jo 

4E^^[^"(^.Mm _,^,(^,^,) - 

Thus we define the Ekman pumping term at the top of the fluid by 
Going back to Lemma [STZt we deduce that 
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3.2. The envelope equation. Now that the term w'"* is defined, there remains to construct w and 
(5it'"* such that m'"* is an approximate solution of equation (|1.2p . We recall that i5m'"*,w'"* are strongly 
oscillating terms, small in norm. Consequently, setting M'"*(t,T) = C{T)w{t), we have 



dtu'"^ + u"* • Vu"* + -63 A It'"* - Ahu'"^ - i^dW"^ 

e 



C 



dtw + u"* • Vw'"' - A^u'"* + - rS^Ju"' + LSu'" 



9tw + Q ( -, w, w ) - Ahw 



where 



+ i [£ (r) dr {C i-T) Su^^Ht, t))] i^^i + S (^t, , 

Q(t, w, u>) = C{-t)F [\/{£{t)w ® /:(r)?i;)] 



and S is defined by 



S(i,r) :=- 



Thus it is natural to choose w and 6u™^ such that for all t, t, 

9tu; + g(T, w, w) - Ahw + C (-r) PS (t, r) + -dr [C (-r) <5M'"*(t, t)] 



(3.12) 



(3.13) 



The quantity £(— t)PI](<, r) has already been computed in Lemma \3l] (see l|3.8p ). Since w does not 
depend on t, the first idea is to average the above equation on a time interval [0, 6], and to pass to the 
limit as 6* ^ oo in order to derive an equation for w. Assuming that the term Ju'"* is bounded uniformly 
in r, we have 



lim / -dr \C (-r) Su'"\t,T)] dr = 0. 



On the other hand, we have already proved in Lemma \3A\ that 



lim ■ 



C i-T) PS (t, r) dr = -5[CB,3, CT,3] 



5bM + j./?4(a) in L^^{[0,oo)t,L\T x E)). 



Moreover, we have 



Q(r, w,w)= 6'(-^'-^'"+^'')" {Ni,w) (Nrn^w) (Nk, (Ni ■ V)7V„,) Nk, 

and it is proved in 0| that if w is sufficiently smooth, 

1 

- y Q{t,w,w) Q{w,w) 
in the distributional sense, where Q is defined by l|1.6p . Hence, for all S > 0, we define as the unique 



(3.14) 



solution in L°°{E, C{[0, oo), H (1 H°^^)) n L'^iE, Ll^{[0,oo), H^^°)) of the equation 



uoEHn H' 



0,1 
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We refer to Proposition 6.5 p. 145 in [4] and to the comments following Proposition ll.2l in the Introduction 
of this paper for existence and uniqueness results about equation (|3.14p . Notice that if cr G L°°([0, T) x 
[0,oo)^ X E,L'^{T^)) only has a finite number of horizontal modes, then S^{a) € L°°{[0,T] x E,H"^^). 
Moreover, the fact that di{Ak) > in the definition of Sb impHes that the Ekman pumping due to the 
Dirichlet condition at z = induces a damping term in the envelope equation. 

• The idea is then to pass to the Hmit in S^{a) as i5 — > when a satisfies (H1)-(H2), using (|1.9p . 
Let us admit for the time being that the last property of Proposition 11.11 holds . i.e. 

3r] > 0, VA e [-1 - 77, -1 + r;] U [1 - ?7, 1 + ??], £x{a) = 0. (3.15) 

Property l|3.15p entails that the sum in the right-hand side of l|3.1ip bears only on the triplets (/ci, ^2,^3) 
such that 



|Afc-l|>r;, |Afe + l|>77, 

which entails 



Consequently, since a only has a finite number of horizontal modes, we deduce that the sum in the 
definition of S^{<7) is finite. Hence S^{<7) converges as 5 —> in L°°{[0, 00) x E; L^(T^ x [0, a))) towards 



kh^Q 

The same property holds when a has an infinite number of horizontal Fourier modes, provided a is 
sufficiently smooth with respect to the horizontal variable Xh and satisfies (H1)-(H2). 

Thus for all To > 0, the source term S^{a) remains bounded in L°°((0,To) x E,H°^^) as (5 -> 0; 
whence is bounded, uniformly in S, in L°°{E,C{[0,To],nn H^^'^) n L'^{[0,To], H^^°)). Moreover, let w 
be the unique solution in L°° {E, C {[0, 00), H n H°'^)) n L°°{E, Ll^{[0, 00), H^'°)) of 



dtw + Q{w, w) - Ahw + ^ '^Sb{w) + vI3St{(j) = 0, 
A standard energy estimate leads to the following error bound, for all T > 0, 

\\w - Ws\\l'=^([o^t]xE.L^) + ||Vh(w - W5)||l~(£;,L2([0,T]xT)) 

< Cu^WSria) - 5^(a)||ioo(s,i2([o,T]xT)). (3.18) 

Thus, when constructing the approximate solution in the next section, we will use the function , 
but we will keep in mind that converges towards w as, 5 vanishes. 

• Let us now turn to the proof of property p.l5p (which is the same as l|1.9p ). Using (H2), we 
choose 770 > such that 

[-1 - 770, -1 + 7/0] CV-, [1 - 770, 1 + 770] C V+. 
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For A G R arbitrary, and for 6* > 0, we have 



L°°([0,T]xT2x_E) 



L°°([0,T]xT2x£;) 



< ||cr — crQ||L=°((0,e)x[0,T]x_ExT2) 
1 



I f e-°'^f'^+"'''~'^^Taa{n)dndT 
Jo Jm. 



L°°{lO,T]xT'^xE) 



< ||cr — crQ||L°°((0,6i)x[0,T]x£;xT2) 



J^aa(Ai) dfi 



i{li-X)e 

< ||cr — crQ||^oo((o,6l)x[0,T]x£;xT2) 



L~(T2) 



+ sup \\^a(T{n)\\^^nQj,,^E^j2) {\V+\ + \V-\) 



R\(V-UV+) 



||jFQCT(/i)||^^QQ ,p]xExT2) djidr. 



Let us now evaluate the last integral when A is close to ±1, say for instance 

|A-1|<^. 
I 1-2 

Then if ^ e M \ {V- U V+), we have l/x - 1| > 770, and thus 

^70 



Im-a|> 



In particular, 



^i{fj.-X)9 _ 2 



< 



\^^-X\9 



< 



c 



Hence, for all 9 > 0, for A such that |A ± 1| < f]o/2, the following inequality holds for all a > 
1 /"^ 

aiT,uj)e-'^^ dT 

° L°^{[a,T]xT^,L^(E)) 

C 

< Ik - O-a||L~([0,9]x[0,T]x£;xT2) + SUp ||^aCr(Ai) II L~ ([O.T] x BxT^.Li (K;^)) 
+ sup ||^aO-(M)llLoo(ro.Tlx_ExT2) (l^+l + l^-l) ■ 

In the above inequality, we first take 9 large enough, so that the left-hand side is close to ||(t(A)||, and 
C sup„ II J^a,(T||/6' is small. Then we let a go to zero, with 9 fixed; we deduce that 

£x[a] = VA such that d(A, ±1) < y . 

3.3. Definition of Su'"*. Once w (or w^) and f'"' are defined, there only remains to obtain an 
equation on (5m'"*. As stated before, Ju'"* is chosen so that equality l|3.13p holds for all r > 0. According 
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to the above computations, this amounts to taking (5m'"* such that 
d 



dr 



£(-t)(5u'"*(t) 



(5m'"' (t) 



(3.19) 



eQ{w, w) - eQ{T, w, w) + S[cb,3, ct.s] - eC{-T)FT,{t, r), 
[Q{w,w) — Q{s,w,w)\ ds 

+ I [S[cB.z.cT.z] - eC{-s)¥^{t,s)\ ds 
Jo 

e£(r) / [Q{w,w) — Q{s,w,w)j ds 
Jo 

+C{t) f [S[cB,3,CT,3]-eCi^s)Fi:it,s)] ds. 
Jo 

Equivalently, Su™^ satifies the equation 

a^^u'"' + L(5u"* = ££(t) [Q{w, w) - eQ{T, w, w)] + C{t)S[cb,3,ct.,3] - eFS(t, t). 

We now derive a bound on the coefficients of (5m'"* : 

Lemma 3.2. LetT>0, N > 0, and let w € L°°{E,C{[0,T],n)) such that 

{Nk,w{t))^0 Vfc,|fc| >7V, Vte [0,T]. 

Let E be given by p.l2p . and 5u''"* by l|3.19p . Then for all k e Z'^, for all 77 > 0, there exists a 
constant C^^k such that for all r > 0, for all €,1^, (3 > such that v — 0{e) and y/evP — C'(l) 

||(7V,,(5M»*(t,r))||^„(j^^^j^^,(^j^ < {e + V^){C,^k+VT). 

Remark 3.2. The above Lemma is stated with a function w having only a finite number of Fourier 
modes, which is not the case for the solution of p.l4p in general. However, when constructing the 
approximate solution in the next section, we will consider regularizations of the solution w of the envelope 
equation l|1.10p . so that this issue is in fact unimportant. 

Proof. We begin with the derivation of a bound for the term 

[Q{w,w) — Q{s,w,w)\ ds 



E 



k.l,m 



Notice that the set {l,m) e Z'' x such that {Nm,w) {Ni,w) ^ is finite, and included in Bn x Bn. 
Moreover, if {l,m) G -Bjv x B^ and A/ + A„i ^ Xk, then there exists a constant aN,k > such that 

I A; + A„i - Afel > aN,k. 

As a consequence, we have 

Nk, / [Qiwit), w{t)) ~ Q{s, wit),w{t))] ds 



< 



OtN,k 



L°'((0,T)xT2x[0,a]x£;)- 



In a similar way, we now derive a bound on the second term in (|3.19p . According to Lemma [XT] we 
have, for all A; G Z^, 

{Nk,C{~s)f^{t,s)) ds^-{Nk,S[cB.3.CT.3]) 
as r ^ c», in L°°([0, T], ^^(i;)). Let r^,fe > such that if r > T,,,fc, then 



- f {Nk,Ci~s)FJ:{t,s)) ds~-{Nk,S[cB,3,CT,3]) 

T Jo ^ 



< 



L^{[0,T],L^{E)) 
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Now, for T < T,i^k, we have 




-S[CB,3,CT,3]- C{-s)PJ:{t,s) ds 



L~([0,T],L2(£;)) 



Nk, -S[CB.3,CT,3] 

' ■ L~([0,T],L2(E)) 
II {Nk, s)) II L=°{[0,T],L^{E)) 

Gathering all the estimates, we infer the inequality announced in Lemma [331 □ 

4. Proof of convergence. This section is devoted to the proof of the convergence result in Theorem 
11.31 In the previous sections, we have already defined boundary layer terms and interior terms at the 
main order. Unfortunately, the sum of those first order terms is not a sufficiently good approximation of 
u^'" . Hence the first step of the proof is to define additional correctors, and thus to build an adequate 
approximate solution. We then derive some technical estimates on the various terms of the approximate 
solution, and eventually we prove the convergence thanks to an energy estimate. 

4.1. Building an approximate solution. The approximate solution is obtained as the sum of 
some interior terms and some boundary layer terms; although we have to construct several correctors 
in order to obtain a good approximation of the function u^^", we emphasize that all terms vanish in 
norm, except the solution of the approximated envelope equation p.l4p . In this paragraph, we build 
the correctors step by step, using the general constructions of the previous sections. At each step, we 
will give some bounds on the corresponding term; these bounds will be proved in the next paragraph. 

• First step. The interior term at the main order. 

We have seen that the interior term at main order is given as the solution of the envelope equation 
(|3.14p . and that when the parameter S vanishes, the envelope equation becomes p.l7p . However, we 
are not able to construct the boundary layer terms at the top for S — 0, and thus we must keep the 
approximated solution of the envelope equation, namely . Moreover, when constructing the corrector 
terms u^^ , i5m'"* , w'"* , we will need some high regularity estimates in space and time on , which are in 
general not available for or w. Thus we introduce a regularization of with respect to the time 
variable, and we truncate the large frequencies in w*. 

Let X € I?(M) be a cut-off function such that 



x{t) = o yte[o,^), x{t) = o vt€(- 



-1], 



x{t) > yte 



For n e N*, set Xn ■— n ^x('/"-): ^ind define, for n, > 0, 



where Pat stands for the projection onto the vector space generated by Nk for |A:| < N. The convolution 
in time is well-defined thanks to the assumptions on the support of %. We have clearly 



lim supljui 
lim sup 1 1 It; 



f'n,JvllL=°([0,T]x_E,L2) 



= 0, 



"n,Ar||L=°(-B,L2([0,T],ffi.0)) = 0. 



n.N 



which 



Moreover, the following result holds, and will be proved in the next paragraph: 

Lemma 4.1. The function wf^ ^ is an approximate solution of l|3.14p . with an error term : 

vanishes in L^([0, T], as n,N —+ oo, uniformly in S. 

Hence we work with wf^ jy instead of w from now on; for all k,s > 0, there exists a constant Cn,N{k, s) 

such that 



ll"'n,JvllL°°(-B,H"=.°°([0,T],_ff»(T)) < Cn,N{k,s). 
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In the sequel, we denote by Cn,N all constants depending on n and TV (and possibly T,uq and cr), but 
not on S. 

• Second step. The boundary layer terms at the first order. 

The boundary layer terms at main order, and Uj^^, are defined in Section O where the function 
w is replaced by wf^ ^. Thus depends in fact on the parameters n,N and S, and u^^ depends on 6. 
Using the results of Proposition 12.21 and the previous step, the following estimates can be proved: 

Lemma 4.2. We recall that v 0{e) and l3y/ei^ = 0{l). Setting 

u^^{t, Xh,z) := u^^{t, Xh,z) + u^^{t, Xh, z) 



ub t, -,Xh, —i= 
£ \/£V 



,BL,ret 



{t,Xh,z) +UT [t, -,Xh, 



we have 



u , zdzU , {z — a)dzU 



BL\ 



lL°=([0,T]xT2x[0,a]x_B) 



N, 



(4.1) 



|L,BL|| ^ ^ ,1/4 

P \\l°^([0,T]xE,H^'''}- "^'^ ' 

Moreover, u^^ is an approximate solution of the linear part of equation i|1.2p . vnth an error term 
hounded in L°°([0,r] x E,L^{T)) by 



The above Lemma follows immediately from Lemma [2711 Proposition I2.2[ Lemma [2731 and Corollary 



Third step. The interior corrector terms u'"* and Su*"K 
We now define the correctors t;'"* and as in l|3.2p - p.3p and p.l9p respectively, taking w 



in (|3.19p . Notice that the boundary conditions cb,3 and ct,3 are of order one in L°° . More precisely, 
using the fact that wf^ has a finite number of Fourier modes on the one hand, and (H1)-(H2) on the 
other, we deduce that 



ll«'"*l|L~([0,T]x[0,oo)xT2x[0,a]) < C (Vi^H W^,Af II L°°([0,T],_f/3) + I^£/3) < Cn,N^V£\ 

moreover, according to Lemma [3?2l 



Vry > 0, Vfc e 3a,,fc > 0, 
Thus we set, for > arbitrary. 



< 77 + C,,^fc£. 



L°°([0,T],L2(£;)) 



\k\<K 

According to the above convergence result, for all if e N, we have 



as e, — > 0. 



L°=([0,T],L2(£;jyl,oc(Y))) 



Moreover, there exists a constant Cn,N,K such that 



''K 



< a 



n,N,K. 



L=°([0,T]x_B,H'i.°°(T)) 



In the rest of the paper, we set 



(4.2) 



3Q 



the following lemma holds: 

Lemma 4.3. Let ^ be the remainder term in the equation on wf^ ^ (see Lemma \4-l\ )- Then the 
function u"** satisfies 

dtu'-^* + -63 A u™* + u™* • Vm^"* - A^u^"* - + Vp = £ ( - ) w + w^r™ + ^"2"" + w's"", 

w/iere wf™ = o(l) m L2([0,r] x x x [0,a]), twf" = o(l) in L^{[0,T] x E,H-^'°), and 
\/n,N, lim sup ||wr'"IU~(£;,-L2([o,T]xT2x[o.Q]) = 0. 

Moreover, 

4to = "0 + 0(1) mL°°(£;,L2(xT)), 
anrf i/iere exists a constant Cn,N,K such that 

II'"™*IIl°°([0,T]xS,M^1'°=(T)) < Cn,N,K- (4.3) 

In the above Lemma and in the rest of the paper, the o(l) means that for all n, N, K, the corre- 
sponding expression vanishes as e, ^ 0, uniformly in S. 

• Fourth step. The boundary layer term at the second order. 

At this stage, we have exhibited a function u'"* (resp. u^^) which is an approximate solution of 
the evolution equation l|1.2p (resp. of its linear part); moreover, the boundary layer term u^^ and the 
corrector i;'"* have been built so that the boundary conditions are satisfied at the leading order. Precisely, 
we have 

d. {<i^.ut) + = Mt, tie) + -i=ac«B,Mc=^ ve) + 

The terms U'y\q^^^^ dcUB,h\c,=^^) '^b,?,\c,=^^ a-nd u^^^*'** are exponentially small, thus satisfy the 
assumptions of Lemma A. 2 in the Appendix; they will be taken care of at the very last step. But in 
general, setting CB,h '■= ^h|z=o~'~'^"K\|z=o' quantity e~^CB,h does not vanish. Thus, we define another 
boundary layer term in order to restore the Dirichlet boundary condition at z = 0. We now have to make 
precise which parts are almost periodic or random stationary in CB,h{t,T). We have 

<U = = -^V.A^i(cT,3) + ^VhA-\cb,3). 

The first term in the right-hand side is clearly random and stationary, whereas the second one is almost 
periodic. Concerning the term (Jm'^*, the situation is not so clear. Using (|3.19p . we write 

<5u'^*(i,r)= e-''^''-Sbk{t,T)Nk, 

\k\<K 



where 



5bk{t,T) := e (^Nk, (Q(w^^^, wf^^^) - Q(s, wf, i^n,jv)) ds^ 
+ (Nk, f S[cB,3, CT,3] - eC{-s)PJ:{t, s)) . 
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According to Lemma 



sup 

te[o,T] 



and 



sup 

te[o,T] 



dt 



Sbk t, 



t 



= 0(1), 



.N- 



Thus we forget the fact that Sbk depends on the microscopic time variable r, and we merely treat (5u^* 
as an almost periodic function. Hence we use the construction of section [2] (see in particular Remark 
2.21 for the random stationary part), and we denote by 6u^^ the boundary layer term thus obtained. By 
definition, 



and 



dtSu^^ + -€3 A 6u^^ ~ vdHu^^ = o(l) in L°^([0,T] x x T). 



Using the same kind of estimates as in Lemma 12. 3^ we deduce that 



• Fifth step. The "stopping" corrector. 

Let us now examine the remaining boundary conditions. 
t> Horizontal component at z = 0: this term is the simplest of all. We have 



and thus, using the same arguments as in Proposition 12.21 we prove that there exists a constant 
C such that 



||fe,/i(i)llH3(T2) < C exp 



C 

\dtSBji{t)\\H3(T^) < — exp 



Since exp {—aj^/ev) = o(l) for all fc G N*, Sb^k satisfies the conditions of Lemma A. 2 in the 
Appendix. 

[> Vertical component at z = 0: we compute 

SbM «(i) +^?^(i) +^"f^(0)|,.o ="T,3|c=^ (*,*/£) +'^^?|^=o(0- 

It is easily proved that MT,3|c=a/VH7(^' V^) satisfies the hypotheses of Lemma A. 2, provided a is 
sufficiently smooth. Concerning (Suf^, we have, according to the assumptions on cr, 

ll^"3|2=olU°=([o,T],L2(£;,//3(T2)) < o(Vei^) + Cn,N ,K{vef^'^ ii < o{e), 

\\dtSu^l"^^o\\L--{[Q,T].L^iE,H3{T^))) = o(l). 

Thus Sb,3 satisfies the conditions of Lemma A. 2. 
> Horizontal component at z — a: 



St At) = + nf^it) + Suf^{t))^^^^ ~ -a (^<, - 



dcUB.HK=^ {t, tie) + d.v^^^;Z{t) + dMl\.=M 
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For all s > 0, we have 



do BL.res 



BL,res 
h\z—a 



L°=([0,T]xS,H»(T2)) 



1 



.N- 



exp 
exp 



0(e). 



(Remember that y = 0(e)-) Thus all terms of the right-hand side are exponentially small as e 
vanishes, and St.h satisfies the conditions of Lemma A. 2. 
[> Vertical component &t z = a: let 

Once again, 5t,3 is exponentially small in all H'^ norms, and thus matches the conditions of 
Lemma A. 2. 

We thus define u^^°^ , given by Lemma A. 2, so that 



stop 



3|z=0 



^B.3, 



^ h\z—a 

stop 

-*3|z=a ~ 



-<5t,3, 



and such that u^*°p is an approximate solution of the linear part of equation l|1.2p , with an error term 
which is o(l) in . Notice that the corrector m^*°p itself is o(e) in . 



We now define 



,app 



,BL 



5u 



BL 



stop 



(4.4) 

= u"' + u™™. (4.5) 

By construction, the remainder u"^™ is o(l) in L°°([0, T], L^(£^ x T) and u^pp satisfies conditions (|1.3p . 
The goal of the next paragraph is to prove that u^pp is an approximate solution of l|1.2p , which allows 
us to conclude in paragraph 14.31 that u'^'" — m^pp vanishes thanks to an energy estimate. 

4.2. Estimates on the approximate solution. We start by proving the lemmas stated in the 
previous paragraph. 

• Proof of Lemma \4Tl\ (Estimates on wf^ j^). 
Remembering (|3.14p , it is easily checked that wf^ ^ satisfies 



Sb{wIj^) + I//3Pjv4('^ * Xn) = 0. 



Thus If ^ jv is an approximate solution of l|3.14p . with an error term jy equal to 

= [(P - Vn)Q{w^,w^)] * Xn + [Qi^NW^, Pjvw^) - Qiw^w^)] * Xn (4.6) 

+ [Q{wlN,wi,N) - Q{Pnw',¥nw') * Xn] + J^/3Pjv4[^ ~ <T * Xn] + i^/3(P " Pn)SU'J)- 

In order to evaluate ^, we need continuity estimates on the quadratic term Q. We recall that Q is 
bilinear continuous from 

L°°([0,r],i/°'i) X L2([o,T],i/i^O) into L^qq, T], if-L"). 

(see Proposition 6.6 in for a proof of this non trivial fact). Moreover, for a,b G n 7i, it can be 
proved, using the methods of [4], that there exists a constant C > such that 



|!0(a,6)lk-i,o <q|a||^//||a||^/f,„||fo||^/'| 



1 1/2 



+C||93a||LHI&llii'l|fe||ffi^ + C||a36|UHl«llL/lla|i;;r,o 
< C (||a||ffi.o||6||ffo.i + ||a||HO,i||5||ffi.o) . 



1/2|, ||l/2 



(4.7) 
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It is easily deduced from the above inequality that the two terms in line (|4.6p converge towards zero; on 
the other hand, the regularity of a entails that S^[(j — a * Xn] vahishes in as rt ^ oo, uniformly in 5, 
together with (P — ¥is[)S^{a). We thus focus on the last term in the expression of ^, which we write 

= [ Qiwi^Nit),PNw\u))Xnit-u)du- [ Qi¥Nw\u),PNw\u))Xn{t-u)du 

Jr Jr 

Q{wi j^it) - ¥Nw\u),¥Nw\u))xnit - u) du, 



and thus, using inequality (|4.7p together with the L°° {[0,T], H^'^) bound on , we infer 

||Q«jvW,<jvW) - Qiw',w') * Xn{t)\\^_,,o 

<C I \\wi^N{i)-^NW^{u)\\ \\Vi,w\u)\\^,^, Xn{t-U) du 

JR 

+ C \\wi^pf{t) -¥NW^(u)\\^^ a\\¥NW^{u)\\jjo,l Xn{t - u) du. 



Eventually, we get 

<C sup \\w'-THw'\\^^ +C sup -Thw'W^^ 

\h\<^ \h\<^ 

where ThW : {t,x) i-^ wit + h,x). The right-hand side of the above inequality vanishes as n — !■ cx3, 
uniformly in 5. 

Thus ^ vanishes as n, — > oo in L^([0,T] x E,H^^''^), uniformly in 6. □ 

Hence we have proved that ^ is an approximate solution of l|3.14p . We now tackle the bounds 
on u'"*. 

• Proof of Lemma lJ^ (Estimates on w™*/ First of all, the estimate (|4.3p is easily deduced from the 
previous bounds on ^, f'"* and Su™^. Thus the main point is to check that the assertions on w. 
i = 1,2, 3, hold true. 

We begin wit 
precisely, we have 

Remembering the definition of Su™^ (see P-lOp l. we infer 

l|w'™™IU~(S,i2([0,T]xT2x[0,a]) 

< II(p-Pk) [q«^,<^)]||^^(^_^,((„,,,^^)) 

+ II(P-Pk) [Q(-,<^,<^)]|L..([0,oo).xi.,L^(MxT)) 

+ ^\\{F-Fk) [^[cB,3,CT,3]]|L^(^_i.([o,T]xT)) 
+ ||(P - IPif)S|lLoc([o,oo)x£;,L2([o,T]xT)) ' 

1{ V = 0{e), and y/ueP — 0{1), all terms vanish as X ^ oo uniformly in e, v, 5. Thus the condition on 
vf^™ is satisfied. 

On the other hand, we have defined w'"* and 5u™^ so that u'"* is an approximate solution of equation 
(fOj) . with an error term which we now evaluate in L^([0,T] x T x i;) + L'^{[G,T] x E,H-^'°). Apart 



i 1 

We begin with the term Wg"™, which is due to the truncation of the large frequencies in 6u™^\ 
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from the one mentioned above, which is due to the truncation of the large spatial frequencies in fe'"*, 
the error term is equal to 



^ (J) <N{t) + m - A, - ^.af )(<Su|* + z;-*)] (^t, 



^ ( - ) ^l^W- 



Let wf^it) := -Ahv"'^ {t,t/e) . Then wf^' is bounded in L'^{[Q,T] x E,H-^-°) by 

V^I|cB,3||L~([0,T]x[0,oo)^x_E,_f/i(T2)) + El'/?!! CT,3 1| ([0,T] x [0,oo)^ x_E,_f/i (T^)) ^ °(-'-)- 

Keeping aside C{t/e) rf^ ^(i), the remaining error terms are bounded in L^{[0, T] x x [0, a] x E) by 

Pt^^^'5*llLoo([o,T]tx[o,f ]^,L2(_ExT)) 
+ ||9tu'"'||L^([0,oo)^,L2([0,T]xT2x[0,a]x£;) 
+ \\^'^K^\\L^{[Q,T]tx[0,^]^,L2{E,H^)) 
+ + ■'^'"*llL~([0.T]tx[0,2:]^,L2(s,Hi)) 

+ ll"'"iL~(B,L2([0,T],Hi))||'5u'^* +w'"*|lL2(£;x-([0.T]tx[0,f]^xT)) 
= 0(1). 

Above, we have used the fact that wf^jy whence v™*, 6v}]^, are smooth with respect to the time 
variable t. □ 

• At this stage, we know that m'"* is an approximate solution of (|1.2p . and that u'^™ is an approximate 
solution of the linear part of l|1.2p . such that additionally u™™ = There remains to prove that the 

function u^pp = m'"* + u™™ is an approximate solution of equation l|1.2p . The core of the proof lies in 
the following Lemma: 

Lemma 4.4 (Non linear estimate on the remainder term). For all n,N, as e, — > with v — 0{e) 
and Py/ev = 0(1), we have 

bupp vu -ru vu -tu vu ||L2([o,rixT2x[o,a]x£;) ^• 
Proof. First, we have 

11/ rem . y\ intll 

11*- llL2([o,T]xT2x[0,a]x£;) 



< 



L2([0,T]xT2x[0,a]x_E) IT II L~ ([0,T] x B, WL^) 

The right-hand side vanishes thanks to the estimates of the previous paragraph. 
The other terms are slightly more complicated. We write 

yint . y^rem ^ ^rem . y^rem ^ ^app . y^rem 

= w'^PP • Vu'^^P + U^PP • V {u^^ + Su^^) 
The first term in the right-hand side is bounded in L^([0,T] x x T) by 

||u^PP|U»||u«*°P|U.([O^T]xiJ,Hi) < C^,N,Ke. 

We thus focus on the second term, which we further split into 
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We have 



Ih/^'PP • V,, (v^^ + (^7/^'"') II 

yh[U -tOU j||i2([o,T]x£;xT) 

< ||u^^''||loo([o,T]x£;xT)I|m^^ + Su^^\\i^2(^[QT]xE,H^.O) 



We split the other term as follows 



T2 Jo 



a/2 



T2 Jo 



a/2 



T2 Jo/2 



T2 Ja/2 



For z > a/2, t > 0, we have 



N 



and thus 



{ev) ^ exp 



/ i^t 



exp 



ca 



JT2 Ja/2 



<C„. 



AT 



(ei/) ^ exp 



exp 



ca 



Similarly, 



a/2 



Jt2 Jo 

2 



< C„,Ar/3'' exp 



< C„,Ar(ei^) ^ exp 



We now evaluate the two remaining terms. The idea is the following: since u^^^ vanishes at the boundary, 
we have 

4PP(z)wCz forz = o(l), 
andu3PP(z)«C(z-a) forz-a = o(l), 

and zdzUB, {z — a)dzUT are evaluated in (|4.ip . Moreover, we can split u^pp into 



L 



By definition of w'"* and m''*°p, the vertical components of each of the three terms in brackets vanish at 
z = and z = a; additionally, the first term is bounded in L°°([0, T] x E, 1^^'°°) by a constant Cn,N,K, 
while the (vertical components of the) second and third ones are respectively of order 



a 



n,jv(x/i^+M'/") and 0(^^/4) +o(e) 
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in L°°{[0,T] X E,H^''^). Once again, the formulation u***°p — o(e) must be understood as 



yn,N,K, lim sup£~^||w'*°P|| = 0. 



As a consequence, we have 







a/2 




Jo 


< 






+ 







\zd, {uB + Su'''^) (t) 



lL2 



\zd, {uB + Su^^) (t) 



Using Hardy's inequality together with the divergence-free property, we infer that 

r ra/2 

JT2 Jo 



lul^^imiuB + Su^'^) it) 

't 



+ Cn.N.K 



t 



L2 



<0{l). 



t 



The term 



T2 J a/2 



wr{t)d,UT{t)Y 



is treated in a similar way. Gathering all the terms, we deduce the convergence result stated in Lemma 

1131 □ 

In the rest of this section, following the notations introduced in Lemma [431 we denote by w™™ any 
term which satisfies 



\/n,N,K, limsup||u;J<""|1^2([o,T]xi5xT2x[o,a]) = 0' 

by any term which satisfies 

Wn,N,K, limsup||w;5<""||^2([o^jx£_^„i.o) =0, 

and by w™™ any term which satisfies 



' s>o 



Vn,iV, lim sup \\wl™\\L^(io,oo)xE,L2ao,T]xT^xlO,a]) = 

According to Lemmas 14.21 14.31 and 14.4} u^pp satisfies an equation of the type 



atw'^PP + u'^PP • Vu^PP + -63 A u'^PP - A,iU^PP - Z/9^uapp 

e 



(4.6 



(4.9) 



(4.10) 



Vp + /:( -^rf^^^+wf^ + uir^ + wS'^^ + oJ^-^^ , (4.11) 
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We recall that the remainder ^, which was defined by l|4.6p . satisfies 



lim sup ||r„^||i2([o,T]xis,if-i.o) = 0. 
Equation l|4.1ip is supplemented with the boundary conditions (|1.3p and the initial condition 



where 6wq and 5wq are such that 



lim sup ||(5wo||ioo(£^2(x)) = 0, 



and Vn,iV, lim sup ||(5woI|l~(e.l2(t)) = 0. 

In order to avoid too heavy notation, we will simply write 

"14=0 = w'o + 0(1). 

4.3. Energy estimate. We now evaluate the difi'erence between u^'"^ and u^pp thanks to an energy 
estimate. The function m'^'" — u^pp is a solution of 



t_ 



Taking the scalar product of the above equation by u^''' — u^pp and using the Cauchy-Schwarz 
inequality, we deduce that for alH > 0, for almost every uj £ E, 

ll\\u^-'^it,u;) - «^PP(t,c.)||i. + \\\u^-''{i.uj) - «-PP(i,c.)||l,.,o 
< / \{{{u^'''{t,uj) - w^PP(t,w)) • V)M^PP(i,w)] • t^) - u^PP(t,u;))| 

+ IK"Xt,^)lli^(T) + IK'"(i,c.)llH-.o + IK-(i,^)|li.(T) 

In the above inequality, we have dropped the term v\\dz{u'^'" — u^PP)|j|2 in the left-hand side. We now 
evaluate the term 



\{{u^^" - M^PP) • V) w^PP • (v!"'" - w^PP)| . 

T2x[0,a] 

First, let us write 

yapp ^ [yint _^ yStopj ^,^BL _^ ^^BLj ^ 

The function uint + u^top is bounded in L°°([0,T] x VFi'°°(T) by a constant C„,Ar; similarly, V,,(u^i^^ 
bvP"^) is bounded in L°°([0, T] x x T). As a consequence, we have 

r int_^yStopj . (ye,i._yapp)| 



/ 




- M^PP) 


JT2x[0,a] 






1 / ^,1^ 




JT2x[0,a] 










«^PP||i2, 
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There remains to derive a bound for the term 



/ 



T2x[0,o] 



Su 



BLl 



the calculations are quite similar to those of Lemma 14.41 We first split the integral on [0, a] into two 
integrals, one bearing on [0, a/2] and the other on [a/2, a]. The term Uj^ (resp. + SvP^) is expo- 
nentially small on [0, a/2] (resp. on [a/2, a]), and thus we neglect it in the final estimate. Moreover, we 
have for instance 



a/2 



J2 



< 



,BL1 



L2 



appi 



L2(T) 



Eventually, we infer that 



T2x[0,a 



3 "3 
2 



,BL 



,BL1 



(^e,!^ _yapp)| 



< C llu^''' - M''PP||i2(x) + CWu^^'' - m^pp||hi,o Wu"'" - u^PP||i2(x) 
Gathering all the above estimates and integrating on E, we deduce that 



d_ 

dt 



,,app| 



|2 

Ll2(£;xT) 



appi 



II^™'"IIl2(£xT) + \\^^2™\\l'^{E,H-^-°) + II"'3''™I!l2(bxT) 



Using Gronwall's Lemma, we infer that for all t e [0,T], 



< C 

+ c 



ll^r™llL2([o,T]x£;xT) + ll'^2''™llL2([0,T]x£;,ff-i-0) + II^S*"" II L2([o,T] x fix T) 

II ' f ^- 

Ir n,Jv||i2([o,T]x£;,H-i'0) + 



(4.12) 



We are now ready to prove Theorem 11.31 Let us write 
't 



u^'^it) wit) 



.appi 



(0 



where 



> the term u"^'" — u^pp satisfies the energy estimate l|4.12p : 



i> the term u^PP(t) - C i 



(t) is equal to u^' 



' and thus vanishes as e, — > in 



L°°{[0,T],L^{E,H^^°)), uniformly in (5 > 0, and for all n,N,K; 
[> the term wf^ — vanishes as n, — > oo uniformly in (5, e,;^ according to the first step in 
paragraph 14. Ij 
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> the term — w vanishes as 5 ^ 0, uniformly in e, z/, according to l|3.18p . 
Let 77 > be arbitrary. We first take no, Nq large enough so that for all S > 0, 6,1^, (H > 0, 



lknn,A'olli°°([0,T]x£;,ff-i-0) ^ 



i5|i2 II S 



..<5||2 



\\l^{E,L2([0,T],H^-°)) - ^■ 

Thanks to l|4.10p . we now choose K > large enough so that for all e,v^(3,5, 

II^3''™IIl2([o^T]x£;xT) - ^■ 

Remembering properties l|4.8p - l|4.9p . we deduce that there exists £01^0 > such that for all 5, for all 
£ < eQ,v < vq with V < Ce and !i\fev < C, 

ll^r™llL2([0,T]x£'xT) - ^' 
1 1 ^,,rem 1 1 2 



"2 \\L^([a,T]xE 
2 



it) 



< 



L°°([0,T],L2(_E,ffi.o)) 

At this stage, we have, for all 5 > 0, for all e, i^, (3 such that < e < Eq and v = ©(e), ^/veP = 0(1), 

2 



L2(_ExT) 



L^(E.H^-O) 



ds 



< Cl] + C\\w^ - w\\i'=°{[0,T],L^(ExT)) + C'llw'^ - w\\j^2^[Q T]xEM^-0) 

We now let 5 ^ in the right-hand side, and we obtain 



C6^ 

e 



w{t) 



L^{ExT^x[a,a]) 

for e,v small enough. The convergence result is thus proved 



L^{E,H^ 



ds < Cri 



5. Mean behaviour at the limit. This section is devoted to the proof of Proposition [L6l Let us 
recall what the issue is: in general, the source term Sxicr) in (|1.10p is a random function, and thus so is 
w. Hence, our goal is to derive an equation, or a system of equations, on E[w]. We emphasize that such 
a derivation is not always possible, because of the nonlinear term Qiw, w). However, we shall prove that 
the vertical average of Wh, denoted by Wh, is always a deterministic function. Moreover, if the torus is 
nonresonant (see l|1.12p ). then w — 0) solves a linear equation, and thus in this particular case we can 
derive an equation for E[w — w]. 

Our first result is the following: 

Lemma 5.1. Assume that the group transformation (6'T-)^gR is ergodic. Let uq G 7i n , and let w 
be the solution of i|1.10p . Set 



Wh 



Wh- 



Thenwh is the unique solution m C([0, oo), L^(T^))nL^P^([0, oo), iJ^(T^)) of the two-dimensional Navier- 
Stokes equation 



dtWh + Wh ■ ^Wh - AhWh H \=\l -Wh + iyf3E [ST{cr)]h = 0, 

av2 V £ 

Wh\t=0 = lIo'^OJi- 



(5.1) 



In particular, Wh is a deterministic 



m 



Proof. Let us recall that if 



J2 kk)Nk e n, 



then 







Thus we have to project equation (|1.10p onto the horizontal modes, which correspond to fcs = 0. It is 
easily checked that 

Ph {Sb{w)) = SB,h{'Wh) = -^Wh, 

and we recall (see ^Uj and Proposition 6.2 in ^) that there exists a function p e L^(T^) such that for 

aWw eH^nn 

Thus we only have to prove that 

Ph{ST{a))^¥.[STAa)], 
almost surely in E. We use the following fact, of which we postpone the proof: if A G M, then 

Moreover, if A = 0, then 

£a[o']=E[(t] almost surely. (5-3) 
Note also that Afc = if and only if ^3 = 0. Remembering l|3.16p . we deduce from (|5.2p and (|5.3p that 



E 


1 


E 


1 

\k',\' 



- Ph[ST{<y)]- 

Thus the lemma is proved, pending the derivation of (|5.2p and l|5.3p . Concerning i|5.2p . the invariance 
of the probability measure mo with respect to 9r entails that 



E[£x[a]] = E[a] lim i / e"*^^ dr, 



and l|5.2p follows easily. Equality IjS.Sp is a consequence of Birkhoff's ergodic theorem (see [21]) • □ 

The first point in Proposition 1 1 .61 follows easily from the above Lemma (together with Theorem [OJ , 
by simply noticing that the sequence 

exp (^-Jl^ w{t) = ^6-^^'=^ {Nk,w{t)) Nk 

weakly converges in i^([0, T] x T x E) towards 

{Nk,w{t)) Nk = wit) = {wh{t),0). 

Afe=0 
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Remark 5.1. Notice that 



mthPh[ST{^)] - -J ^bE [rot^a] . 
V a 

Hence we recover the result of l&j: the vorticity (j> :— rothWh is a solution of 



dt(f> + Wh ■ Vh(f> - ^h<ii + -<j> = [rothcr] . 

aV2 \ £ \ a 

From now on, we assume that the torus is nonresonant (see (|1.12p ). Consequently, with w — {wh, 0), 
we have 

Q{w — w,w — w) = 0. 
Moreover, using (|5.2p - l|5.3p . it is easily checked that 

Setting u — w — IV, we deduce that u solves a linear equation, namely 

dtu + 2Qiu,w) - Ahu + <^Sb{u) + lypSricr) - i^l3E[ST{a)] = 0. 
Since w is deterministic, we have 

E[Q{u,w)] =Q{E[u],w). 

Hence we can further decompose u into w + u, where w is deterministic and does not depend on a, and u 
is random with zero average. The precise result is stated in the following lemma, from which Proposition 
11.61 follows immediately: 

Lemma 5.2. Assume that the hypotheses of Proposition FTTgl hold. Then 

W — w + w + u 

where: 



the function w is deterministic and satisfies (|5.ip ; 
the function w is deterministic and satisfies 

dtw + 2Q{w, w) - Ahw + ^Sb{w) = 0, 

• the function u is random, with zero average, and satisfies 

dtu + 2Q{w,u) - AhU + ^Sb{u) + lypSria) - i^f3E[STia)] = 0, 
"|t=o = 0. 

Appendix A: convergence of the family ctq. 

Lemma A.l LetT > 0. Assume that a e L°°{[0,T] x E,C{R)) n L°°{[0,T] x x E). Then for all 
T' > 0, 

CTa-cr^O in L°°{{0,T) X {0,T') X E) as a ^ 0. 



Proof. By definition of aa , we have 

o-c((t, r, w) = / e:iq){-a\X\ - a\s\)e^^'-'^^''^a(t,s,uj) ds dX 



27r „ 

1 /■ 1 

— / cxpf—alr + asl) -o{t,T ^ as,u}) ds. 



42 



Consequently, 

If 1 

a{t,T,uj) — aa{t,T,uj) — — / exp(— ajr + asj) :t [a{t,T,Lu) — a{t,T + as,(jj)] ds 

Jr l + s^ 

If 1 

H — a{t,T,LLj) / [1 — exp{—a\T + as\)] ^ ds. 

TT Jr 1 + S 

The convergence result of Lemma A.l follows easily. □ 

Appendix B: proof of Proposition fTTTl Let A g M be arbitrary, and let G L'^{E). 
Consider the probability space 

:=^x [0,27r), P^:=P^^, 

where fi is the standard Lebesgue measure on [0, 27r]. Let us define the following group of transformations, 
acting on {Ex,Px) 

T/(w,(^) (6*^^,^- At mod27r), r e M. 

Then it is easily checked that is measure-preserving for all r e M. And if T > 0, we have, for all 

if e [0,27r], 

Jo 



-tip 



r^{T,\u;,^)) dr, 
Jo 



where the function ^' e L'^{E\) is defined by 

Hence, according to Birkhoff's ergodic theorem (see [21]), there exists a function '^^ G L^{Ex), 
invariant by the group of transformations {T^')_^^^, such that 







T . 

Pa - almost surely in E\ and in L'^{E\). Moreover, the function 

clearly does not depend on (p by construction. Hence, we set 

:= e'"^^^{Lu, if) y{uj, if) e Ex 

and we have proved that 



^ Jo 



almost surely in u and in L'^{E). 

Now, since is invariant by the group (7^^)^^^ and <I>^ does not depend on we have, almost 
surely in w, 

= e-^('^--^^)^'^(6l^w, ip-\T mod27r) 
This completes the proof of Proposition 11.11 
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Appendix C: the stopping Lemma. 

Lemma A. 2 (Stopping condition) Let To > 0, and let 6b, St e L°°{[0,To],H''^{T^)) be two families 
such that 



JiSr. 



3 - 5B,3)dXh = 

and such that as s 0, for A — T, B, 

-Pa||l~([o,t],_h"1(t2)) 0, \\Sa\\l°°{io.t],h^{t^)) and ||9t(5A||L°°([o,T],//i(T2)) ~^ ^■ 
Then there exists a family u^*°p e L°°([0, T], L^(T)) with V • w = such that 

u'TJo = Sb, ul*°l^ = St,3 and d,u'^°l^ = Sr^h 



and such that as e —i- 0, 

-||u**''''||ioo([o_T],L2) 0, 



£ 



-> 0. 

L~([0,T1,L2) 



For a proof of the above Lemma, see [5]. 
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